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Multiple resolution analysis of two dimensional structures composed of randomly adsorbed pene- 
trable rods, for densities below the percolation threshold, has been carried out using box-counting 
functions. It is found that at relevant resolutions, for box-sizes, r, between cutoffs given by the 
average rod length {£) and the average inter-rod distance n, these systems exhibit apparent fractal 
behavior. It is shown that unlike the case of randomly distributed isotropic objects, the upper 
cutoff n is not only a function of the coverage but also depends on the excluded volume, aver- 
aged over the orientational distribution. Moreover, the apparent fractal dimension also depends 
on the orientational distributions of the rods and decreases as it becomes more anisotropic. For 
box sizes smaller than (I) the box counting function is determined by the internal structure of the 
rods, whether simple or itself fractal. Two examples are considered - one of regular rods of one 
dimensional structure and rods which are trimmed into a Cantor set structure which are fractals 
themselves. The models examined are relevant to adsorption of linear molecules and fibers, liquid 
crystals, stress induced fractures and edge imperfections in metal catalysts. We thus obtain a dis- 
tinction between two ranges of length scales: r < (I) where the internal structure of the adsorbed 
objects is probed, and {£) < r < r\ where their distribution is probed, both of which may exhibit 
fractal behavior. This distinction is relevant to the large class of systems which exhibit aggregation 
of a finite density of fractal-like clusters, which includes surface growth in molecular beam epitaxy 
and diffusion-limited-cluster-cluster-aggregation models. 



I. INTRODUCTION 

The structure, packing and phase transitions in media composed of rod-like particles are of great importance in a 
large variety of physical systems Jl]-§|]. Such systems include liquid crystals jfj), slip lines in stressed metals and rocks, 
line-fractures in catalysts |7]-[|], fibers in paper sheets |Io| , pT| ], and rod- like metal particles |l^ , |l3| ]. Unlike systems 
composed of nearly isotropic particles, for rod-like particles the orientational distribution plays an important role 
in addition to the positional distribution jH- Systems of adsorbed particles can be classified to ones of penetrable 
^5|-|l7| vs. impenetrable particles 18 1S|| . In this paper we will focus on the case of penetrable rods. 

Many of the systems of interest in this context are essentially two dimensional and can be studied by a model 
in which M penetrable rods are randomly distributed on a two dimensional surface of size L x L. Each rod is of 
length I and width w. The effective two dimensional coverage is given by 772 = p ■ £w, where p = M / L 2 is the rod 
density. It provides a measure of the fraction of the area which is covered by the rods, where the reduction due to 
overlaps is ignored. Clusters of overlapping rods start to form at low rod density and increase in size as the density 
increases. Systems of this type exhibit a percolation transition at rod density p c , above which there is an cluster 
which spans across the entire system. The percolation transition for systems of rod-like particles has been studied 
extensively in recent years [[l4 16 l^]20| -p4||. A criterion for the percolation density p c was obtained and examined for 



a variety of systems with different distributions of rod sizes and orientations. According to this criterion, percolation 
for isotropically oriented and randomly distributed rods occurs at rod density p c = 1/A exc [[HJ. Here A exc is the 
excluded area for placement of the center of a rod given that it may not intersect a rod already on the surface. This 
area, averaged over the isotropic orientational distribution is JlTj : 
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A exc = 2(l + ^)w£+-(£ 2 +w 2 ). (1.1) 

Note that, since the excluded area A exc , is not simply related to the rod area £ ■ w, the percolation threshold cannot 
be expressed as a function of the coverage 7/2 alone. This distinguishes the rods from isotropic objects such as disks 



for which the excluded area is related to the disk area by a constant factor |17|. This difference is especially important 
in the case of large aspect ratio of the rods. In particular, for w = 0, the coverage vanishes but the excluded area 
does not. 

The properties of the infinite percolation cluster at the transition have been studied extensively for lattice systems 
as well as for continuous systems such the the ones considered here H^]2(|. In particular, its fractal nature was 
established. The fractal dimension, which can be expressed in terms of other critical exponents was found to be 
universal ^,^]. Properties such as correlation functions and cluster size distributions were also examined [^8|-^0|. 
In addition to percolation systems, spatial fractal structures have been observed in a great variety of systems in 
physics |3l|-|36|l and chemistry |57j. 

In this paper we perform multiple resolution analysis of systems composed of randomly adsorbed rods below the 
percolation threshold p c . In many empirical systems of this type, where adsorption and aggregation processes take 
place, fractal structures appear over a range of length scales for densities below p c fl38| . Moreover, a broad spectrum 
of fractal dimensions has been observed experimentally in such systems. To gain a better understanding of empirical 
fractal structures composed of rod-like particles we study these structures using multi-scale analysis, based on the 
box-counting (BC) function. In this analysis we examine the fractal properties of the entire system rather than 
a particular cluster. This procedure is the one most commonly used in experimental studies of adsorption and 
aggregation phenomena. In our analysis, we concentrate on a range of rod densities below the percolation threshold, 
where the area fraction covered by the rods is small. We identify two ranges of length scales, in which the observed 
dimension of the system is different from the space dimension. For length scales larger than the typical rod length, 
the results resemble those obtained for randomly distributed spheres |39|j . An apparent fractal behavior is observed 
within a range of length scales between the rod length and the average distance between adjacent rods, which can span 
up to two decades. However, as was found before in percolation systems, highly anisotropic particles exhibit new and 
important features which do not appear in the study of isotropic particles such as spheres. In particular, we observe 
that, unlike the case of spheres, the upper cutoff (and thus the range over which apparent fractal behavior appears) 
cannot be expressed in terms of the effective coverage alone, but also depends on the excluded area. Moreover, we 
find that the apparent fractal dimension (FD) depends on the orientational distribution and decreases as it becomes 
more anisotropic. For length scales below the rod length, the BC function is dominated by the structure of the single 
rod whether simple or itself fractal. The class of systems which exhibits these two ranges includes the broad class of 
systems in which fractal-like clusters are adsorbed on a surface. To examine the fractal properties of such systems 
we study a simple one dimensional model in which the randomly distributed objects are Cantor sets. The paper 
is organized as follows. The model of randomly distributed rods is introduced in Section II and the box-counting 
function is examined. Rods of zero width are studied in Section III, randomly distributed Cantor sets in Section IV 
and a summary is presented in Section V. 



II. THE GENERAL MODEL OF RANDOMLY ADSORBED RODS 



In this Section we study a rather general model in which M rectangular rods are randomly adsorbed on a square 
surface of area L 2 . They are adsorbed with no correlations and are allowed to overlap with each other. The lengths 
of the rods < £j < £ m ax, i = 1,...,M are independently picked from the distribution Pe(£) and the widths 
< Wi < w max , i = 1, . . . , M are picked from the distribution P w (w) . The orientations < 9i < tt are picked from 
the distribution Pg(9). 

In order to examine the apparent fractal nature of the resulting random structures we apply the BC procedure. In 
this procedure one partitions the surface into a grid of mesh size r ■ L (0 < r < 1) and counts the number of boxes, 
N(r), intersected by at least one rod, as a function of r. An analytical solution for the BC function N(r) will now be 
obtained, and its slope examined in the different domains. 

Imagine the surface of area L 2 to be initially empty and arbitrarily choose one of the grid-boxes of size r ■ L. Denote 
the excluded area for placement of the center of rod i, given that the rod may not intersect the chosen box, by 
S r (£i, Wi, 9i). This area is shown shaded in Fig. |l|(a). Simple geometrical arguments yield: 

S r (£i,Wi,6i) = (r • L) 2 + ii ■ Wl + r ■ L(cos<9 4 + sin^)(^ + w t ). (2.1) 

The two-fold rotation symmetry of the rods indicates that the angles Qi are distributed in the range < 0i < n. 
However, the four-fold symmetry of the box makes the problem invariant under the transformation 9^9 + ir/2. 
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Therefore, we can simplify the calculations by considering only the range < Oi < ir/2 with no effect on the results. 
The angle and side-lengths distributions are thus normalized so that: 
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d£Pe(£) = / dwP w (w) = 1. 
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(2.2) 



We denote quantities averaged with respect to these distributions by angular brackets: (■■•). 

Following Ref. 15 1 we consider the probability q\ for random placement of the first rod without intersecting the 
box. This probability, which is proportional to the free area, is given by: q\ — [L 2 — S r (i\, w\, 9i)]/L 2 . The next 
rod is placed with new random angle and side-lengths, so that the probability for two successful placements is: 
<Z2 = [1 - S r (*i,u;i,0i)/L a ][l - S r (e 2 ,w 2 ,9 2 )/L 2 }. Clearly, q M = nf=i [l - #r(4 w h 8^/L 2 ] is the probability of 
placing M rods without intersection with the chosen box. Thus, the probability of at least one intersection after M 
placements is pm = 1 — qM- Since the total number of boxes is 1/r 2 , for a given realization of angles and side-lengths 
the expected number of intersected boxes is pm/t 2 , or: 
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This expression still has to be averaged over the side-lengths and angle ensembles: 
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where the last equality foll ows since the 6i are independent, identically distributed random variables, and so are li 
and Wi. Finally, using Eq. (2.1) this can be expressed in terms of averages: 
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(2.5) 



This BC function contains the information required for the multi-scale analysis of the system. In order to proceed 
with the scaling analysis, we turn next to the identification of the cutoffs, which separate between length scales in 
which N(r) is dominated by the distribution vs. the internal structure of the rods. These cutoffs are determined 
by the effective dimensions of the rods, when projected onto the grid of boxes. Now, for a given rod i, which is 
oriented at an angle 9i, the average between the projections along the x and y axes (parallel to the box sides) is 
£i(cos9i + sin#i)/2 for the length and Wi(cos6i + sin#i)/2 for the width [Fig. [|(b)]. To obtain the effective cutoff, 
this should be averaged over all angles and side- lengths. We thus expect to find cutoffs at: 
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Unlike the lower cutoffs which arc determined by the rod dimensions, the upper cutoff (ri) is given by the average 
distance between adjacent rod sides [p9|: 
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(2.7) 



To gain insight into the behavior of Eq. ( [2.5D we display it in Fig. |2j for the case of a random distribution of 
identical rods of length I = d\ and width w = d 2 [d\ 3> d 2 ) with isotropic orientations. We chose the case of narrow 
rods, since virtually all of the experimentally relevant cases mentioned above, belong to this category. In the next 
section we return to analyze the non-isotropic and polydispersed cases. The non-trivial apparent fractal region, is 
clearly seen in Fig. 2, between the predicted cutoffs of log 10 (7"o) = —4.2 and log 10 (ri) = —1.5. It is the range which is 
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typical to all (fractal) resolution analyses, and which we claim, reveals apparent fractality. The apparent FD for this 
region (see derivation below) is D = 0.06. In Fig. ^, all four different approximately straight line regions are clearly 
identified, corresponding to the presence of the three cutoffs. The range r < r™ , is usually uninteresting from the 
experimental point of view because it probes the structure of the building block itself . It has a -slope (dimension) 
of 2, namely the dimension of the underlying plane. When < r < Tq, the resolution of observation is coarser, and 
the rods appear as one-dimensional objects, reflected in the -slope (dimension) of approximately 1. The width of this 
range depends on the aspect ratio of the rods and is approximately log 10 ((^)/ (w)) in decades. 

Having the cut-offs at hand, we now perform a scaling analysis in the region between them. This will be done by 
applying the standard fractal procedure 

log 10 iV(r)~-£Mog 10 (r) (2.8) 

where D is the apparent fractal dimension (FD). As stated above, the range which has been the focus of attention is the 
one in between the cutoffs r$ < r < r\. The apparent FD, given by the slope in this range, is a non-universal dimension, 
the magnitude of which depends on the coverage. In order to obtain expressions for the coverage dependence, it is 
convenient to define effective one dimensional (ID) and two dimensional (2D) coverages. The effective 2D coverage is 

, a = M<^. (2.9) 

while the effective ID coverage is 

m=v{ + v? (2.io) 



where iji = vM?-q and rff = (note that as the rod widths approach zero 771 — * rj[). The significance of the 

effective ID coverage is that it provides a unified measure of coverage which is independent of the space dimension. 

The width of the apparent fractal range, can be estimated by A e = log 10 (ri) — log 10 (^o)- I n the case of large aspect 
ratio ((id) < (£)): 

A e = lim (logn - logr£) « log(l/^ - 1). (2.11) 

(w)/(t)-*0 

Since in the limit of large aspect ratio rji cx (p ■ -Ag^c) 1 / 2 up to a factor of order 1, it is found, from the criterion 
of Ref. JTtJ (p c = 1/A exc ) that the width A e is nearly zero in the vicinity of the percolation threshold, and increases 
as the density decreases below this threshold. This estimate for A e may seem to suggest that one can increase the 
scaling range at will by decreasing the coverage. However, in addition to the width of the range between the cutoffs, 
the quality of the linear fit within this range measured by the coefficient of determination R 2 should be considered. 
One can limit the range of linearity by imposing a lower bound on R 2 . Obviously, the range decreases as R 2 increases. 
We thus conclude that the two cutoffs limit the width of the linear range for high coverage while the R 2 criterion 
limits it for low coverage [see Ref. |39) and compare Figs. ||(a) and ||(b)]. The net result is that the range of scaling 
is in fact restricted to 1-2 decades. 

The apparent dimension is found by calculating the logarithmic derivative of N(r) at the estimated middle point 
of the linear range. This point is situated at: r m = yr\ • 7q. Assuming (w) < (I) this yields for the apparent FD, in 
the M — > 00 limit and at constant coverage, the general equation: 
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(2.12) 



This general formula for the FD will be used in the analysis to follow. The effects of changing the coverage on the 
apparent dimension, are shown in the inset of Fig. |^ for dx/d® = 10 3 : the FD increases monotonically with coverage, 
but does not reach 2 for the relatively narrow rods considered. 

The dimension of the rods themselves can also be obtained from the box counting function. In order to do so one 
needs to identify the appropriate length scale and measure the slope of the box counting function N(r) on the log — log 
plot. This slope should be obtained for a length which is the geometrical average of the effective cutoffs associated 
with the average rod length (£) and width (w) according to r = \J ■ r l Q . In the zero- width rod limit dild\ — ► this 
dimension approaches 1. As we are focusing especially on the experimentally relevant case of narrow rods, and as it is 
evident from Eqs. (|J), fl2.10| ) that both the FD and the range A e are only marginally dependent on the rod width, 



we consider in the next cases, for simplicity, rods with zero width. 
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III. APPLICATIONS TO RODS OF ZERO WIDTH 



In the present Section we consider for simplicity the case of rods with zero width. Some comparisons with simulations 
will be presented, as well as applications of the general theory to specific distributions of interest. 



A. Isotropically Oriented Identical Rods 



As a first case for the zero-width rods, let us return to the previous example of identical rods with isotropic 
orientations and impose zero width (w — 0). For the isotropic [Pe(0) = 2/7r], monodispersed \PeU) = 5(1 — d)} case, 

(EL 



the BC function and FD for rods of length d and zero width follow directly from Eqs. (2.5), (|2.12) respectively: 
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In Fig. || we present this analytical result for the rod-BC function along with numerical simulations, for two coverages. 
The agreement between theory and simulations is excellent over the entire range of box-sizes. The ranges of apparent 
fractality are brought in the insets; one would expect the wavy nature of the line at the lower coverage, to be smeared 
out by noise in typical experimental situations. Comparing to the BC function for finite-width rods, which appears 
in Fig. ||, one observes that, as expected, there are now only two cutoffs and correspondingly three (approximately) 
linear regions of slope ~ 1, D < 1 and 2. 



B. Anisotropically Oriented Identical Rods 



Anisotropically oriented elongated particles appear in a wide variety of systems, notably in liquid crystals. In 
order to investigate the effect of anisotropy on the apparent fractal properties we consider here the following angular 
distribution, normalized for < 9 < ir/2: 
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(3.3) 



In the limit n — > this corresponds to a uniformly random distribution, whereas for n — > oo, to perfectly aligned 
rods. T he r ods are assumed to be of equal size d. The BC function can be found by calculating the angular averages 
of Eq. (2.5). Using the identities (n > 0): 
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we find: 

7„ = ((cos I 
This yields for the BC function: 



-sine*)) 
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^ r(l/2 + n) \T{3/2 + n) + V^(l + 2n) 
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In the n — > limit one retrieves the result for uniformly randomly oriented rods , Eq. (3.1), whereas for n — > oo one 
finds N(r) = {1 — [1 — r 2 — rd/L] M }/r 2 . For the cutoffs we have from Eq. (^.6[): 



(3.4) 
(3.5) 

(3.6) 
(3.7) 
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and for the coverage: 

r& = VM ln -±- (3.9) 



The FD for the present case is found by substituting r/f from Eq. 
taking rj2 = 0. The effect on D of changing n at constant number and size of rods is shown in Fig. |4|: the apparent 
FD decreases as the rods become more parallel, i.e., as they "cover space" less effectively. For isotropic objects 
we have shown that the apparent FD depends essentially only on the coverage [^9|. As seen here, for anisotropic 
objects, the FD depends on an additional parameter (the degree of anisotropy). This feature may be relevant to many 
experimental systems which exhibit anisotropic distributions of rod-like particles, such as liquid crystals and paper 
fibers. 



3.E) in the general expression Eq. (2.1 2|) and 



C. Polydispersed Rods 



We will now explore the effects of polydispersivity in the rods size. In particular we will consider a power-law 
distribution of the rod lengths. 



1. Model 



In Ref. [ p9[ a variety of narrow size distributions were examined. It was observed that such polydispersivity does 
not alter the basic observation of an apparent fractal regime between cutoffs. The apparent FD of the correspond- 
ing monodispersed distribution was only slightly modified. An important distribution function found in numerous 
experimental cases is the power- law distribution of sizes fififl . To obtain such a distribution we choose rod lengths 
from an iteratively constructed Cantor set, containing 2™ segments of length 3~ n L in the n th iteration. Assuming 
rod lengths are chosen uniformly from among these segments, the probability of choosing a segment of length 3~ n L 
is P(£/L = 3~ n ) = 2 n /Z, where Z — ~^2k=i 2 fc i ^ m being the maximal iteration number in the construction of the 
Cantor set. The average rod length £ is then given by: 



(£) = f^£ P(£/L = 3- fe ) = i ^ 2 k 3- k L = 1 L 



k=l k=l 

Expressed in terms of lengths one finds for the distribution: 
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(3.11) 



D c being the FD of the Cantor set. Thus the length distribution indeed satisfies a power-law 41 . Note that P(£) 
is a discrete distribution with allowed lengths of I = 3~ n L, and that (£}, for example, would be different if £ could 
assume any value between 3~ km L and L/3. 



2. Scaling Analysis for Power-Law Dispersed Rods 



We assume t hat the rods are uniformly randomly oriented, so that (cos 9 + sin ( 
have from Eq. (pT3): 



4/tt, < 9 < tt/2. We thus 
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with (£} given by Eq. (3.1C). As can be seen in Fig. |], for as few as M = 62 rods, Eq. ( |3.12 ) is in good agreement 



with the simulation results for rods with a power-law distribution of lengths. Thus, we find that the BC function can 
be derived analytically for a power-law distribution of particle dimensions, and depends essentially only on the first 
moment of this distribution. The apparent FD can also be found as before, as the slope between cutoffs determined 



by the average rod length. From Eq. (2.6), one gets in the zero-width rod limit: 



2(t) 



7T L 



The effective ID coverage is now given by: 



ri = —== — . (3.13) 
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The FD is found from Eq. ([3.2]), with the present r/f. The prediction of this formula is compared in Fig. || with 
a linear regression in the range set exactly by the above cutoffs, and is in good agreement. Thus the FD is still 
determined essentially only by the ID coverage. This is a non-trivial result, since a power-law distribution is poorly 
described by its mean, yet this is essentially the only distribution-related quantity needed to express the BC function. 
It can be understood intuitively as follows: since the BC function counts the total number of occupied boxes, it is 
approximately proportional to the total length of all rods. This quantity is well described by the number of rods times 
the average rod-length, namely the ID coverage. In the next Section we consider the effect of endowing the rods with 
an internal (fractal) structure. 



IV. RANDOMLY ADSORBED CANTOR-RODS 



So far we have dealt with random distributions of objects which are not themselves fractals. This was reflected 
in the BC function at resolutions below the lower cutoff, by an integer slope. In this Section we consider a model 
of randomly deposited rods which are all Cantor sets of FD=Z? C . A new feature expected in this case, is that for 
perfect (i.e., not truncated) Cantor sets, a non-integer slope of —D c should appear below the lower cutoff, in contrast 
to the cases considered so far. The main motivation for considering a model of Cantor-rods, however, is that it 
mimics a large class of experimental systems where a set of fractal objects is randomly adsorbed on a surface. For 
example, DLA-like clusters growing simultaneously from several nucleation centers |{42|| and cluster-cluster aggregation 
experiments p3]. In such systems we expect an interplay between the FD of the fractal objects and the apparent FD 
induced by randomness. In particular, it might be difficult to disentangle the respective slopes of the BC function if 
the objects are fractal over a small range. 

For simplicity we will now consider the one dimensional case where M Cantor sets (rods) of total length d each, are 
deposited on a line of length L. As in Sec. Ill, let n denote the number of iterations in the construction of the Cantor 
sets. Thus there are 2™ segments of length 3~™ in each rod. As before, the rods are fully penetrable to each other. In 
order to find the BC function for the resulting set, we focus on a single, arbitrary box of length r ■ L and calculate the 
excluded length for placement of the center of a rod, the condition being that the rod does not overlap with the box. 

Let denote the excluded length for placement of a single n th iteration Cantor-rod, when the box-length satisfies: 
3 < r ■ L < 3~ k d. To find this function, consider first the case n = 0, i.e., the case of gapless rods of length 

d. This is nothing but the ID version of equi-sized rods considered in Sec.||, and the excluded length is clearly 
Oq = r ■ L + d. When n = 1 there is a gap of length d/3 and one must distinguish between the cases r ■ L > d/3 
and r ■ L < d/3. In the former, the resolution of the boxes is insufficient to notice the presence of the gap, namely, 
if the box overlaps with the gap it necessarily touches at least one of the two rod-segments as well. In this case, 
therefore, the excluded length is again Qq = r ■ L + d. However, when r ■ L < d/3 a new situation arises: the box can 
fully overlap with the gap. This is equivalent to having two rods of length d/3, each contributing r ■ L + d/3 to the 

excluded length: fij 1 ' = 2(r • L + d/3). When n = 2 there are gaps of length d/3 and d/9, so that three cases arise: 
(1) r ■ L > d/3, (2) d/9 < r ■ L < d/3, (3) r ■ L < d/9. The first two do not differ from n = 1 since the the boxes 

are not small enough to resolve the d/9 gaps: fl — r ■ L + d, Sl\ — 2(r • L + d/3). The third case is equivalent to 
having four rods of length d/9, each contributing an excluded length of r ■ L + d/9: = 4(r • L + d/9). The general 
case should now be clear: 
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(4.1) 



It is convenient to express the index k satisfying the constraint in Eq. (4.1) as: 

k 



[ log[d/(r-L)] J 



log(3) 



(4.2) 



(where [x\ is the largest integer smaller than x). Now suppose r is given, choose an arbitrary box, and place a 
Cantor-rod at random on the line. The probability q\ that the rod does not intersect the box is the relative available 
length, i.e., q\ = 1 — fij^ /L. For M independently placed Cantor-rods the probability that none intersects the chosen 
box is qM — qf 1 , and the probability of at least one intersection is: pm = 1 — 9m- When multiplied by the total 
number of boxes (1/r), this yields the expected number of intersected boxes: 
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(4.3) 



where k is given by Eq. (4.2). This expression for the BC function is plotted in Fig. || for fifth generation Cantor- 
rods. The remarkable feature in comparison with randomly positioned full rods is that the lower cutoff has shifted 
to the left, now reflecting the size of the smallest segment in the Cantor-rods. At smaller resolutions the boxes again 
"see" ID objects and the slope of the BC function is — 1. At a resolution of r = d/L (rod-size), there is a smooth 
transition into the regime of apparent fractality associated with the random distribution of Cantor-rods. We believe 
that such behavior is typical of experimental situations where a random distribution of limited-range fractal objects 
is observed. However, unlike the present "clean" case, it may be much more difficult to separate the two regions in 
actual experimental data. 

It should further be remarked that as can easily be checked, Eq. (4.3) yields a logarithmic derivative for N(r) of 
— D c in the limits n — > oo and r — > 0. This implies the interesting result that regardless of coverage, a random ensemble 
of Cantor sets has exactly the same FD as a single Cantor set. Fig. ^ shows that in the range between the size of the 
smallest segment and the size of the set, the same conclusion holds for truncated Cantor sets. Clearly, our methods of 
calculation can be extended to other iteratively constructed fractal sets, and the conclusions above regarding the FD 
should therefore remain valid in such cases. We conjecture that they remain valid also for non-iteratively constructed 
fractals. 



V. SUMMARY 



We have performed a multiple resolution analysis using box counting functions to structures composed of randomly 
adsorbed rods. Such structures appear in a large variety of adsorption phenomena and in many physical systems 
including liquid crystals, aggregates of linear molecules and fibers and line fractures. The scaling properties of these 
systems are determined by the particle size (length and width) distribution, orientation distribution, as well as by 
correlations in the positions and orientations. In processes such as paper fiber sedimentation particles tend to aggregate 
into dense regions and create inhomogeneities ]l0| , p[ . In other systems, rods cannot overlap giving rise to a maximal 
jamming density Jj^EJ]. 

We have studied the case where rods can overlap and there are no positional or orientational correlations between 
them. We examined various size distributions and orientation distributions, found an analytical expression for the 
box counting function and compared the analytical results to numerical simulations. 

For rod densities below the percolation threshold, we identified two interesting ranges of length-scales in which the 
box counting analysis gives rise to non trivial scaling properties. In the range of length scales between the typical 
rod length and the typical distance between adjacent rods, the box counting function is determined by the positional 
and orientational distribution of the rods in the plane rather than the structure of the single rod. This gives rise 
to an apparent fractal behavior over a finite range of up to two decades. Unlike the case of randomly distributed 
isotropic objects (disks), the range of length scales over which apparent fractal behavior is observed depends not 
only on effective coverage but also on the excluded area (averaged over the orientational distribution) . Moreover, the 
apparent fractal dimension depends on the orientational distribution and decreases as it becomes more anisotropic. 
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For length scales smaller than the typical rod length, the box counting function is determined by the internal 
structure of the rod. In case of ordinary one dimensional rods this gives rise to a dimension of 1. However, this 
result is more general and in case that the rods are trimmed into Cantor sets the FD which is observed in this 
range is equal to the FD of the single Cantor set. The distinction between these two ranges of length scales applies 
for a very broad class of systems which exhibit nucleation of a finite density of fractal-like clusters. The majority 
of spatial fractals in the physics literature belong to this class [p8[ . A finite density of fractal- like clusters appears 
in models such as diffusion-limited-cluster-cluster-aggregation p3 ] and in experimental systems such as molecular 
beam epitaxy where diffusion-limited-aggregation-like clusters nucleate at a finite density |42j| . In these systems the 
fractal-like structure of the clusters results from complex stochastic dynamics. The cluster density is determined by 
parameters such as the temperature and deposition rate. The typical distribution is not exactly Poissonian due to 
effective repulsion between clusters. However, the distinction, emphasized here, between smaller length scales where 
the fractal properties are dominated by the single cluster and larger length scales where the distribution is dominant 
still applies. We thus predict that experiments involving random distributions of fractal-like objects, will reveal a 
cross-over from object-dominated to distribution-dominated fractal behavior. 
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FIG. 1. (a) Excluded area for the placement of the center of a rod of length i = d\ and width w — d-2 inclined at angle 6 
relative to the horizontal axis. If the center falls outside the shaded area, the rod will not intersect the box of side r; (b) At 
the orientation shown, the projections of the sides, d\ and di, on the horizontal axis are d\ cos 8 and efe sin# respectively (the 
sum of the two projections is shown by the thick line). Th e projections on the vertical axis are d\ sm6 and efe cos# respectively. 
Averaging over the two possibilities and 8 yields Eq. 



FIG. 2. The BC function for rods (Eq. 2.5). Four approximately linear regimes indicate the presence of three cutoffs. As 



can be seen, the locations of these match the prediction of Eq. (2.6) [ logi n (r ) = —4.2, log 10 (r™) = —7.2, and log 10 (ri) = —1.5]. 
The inset shows the apparent fractal dimension D, given by Eq. (2.12), (for the range of length scales between r and n) 



plotted as a function of the effective ID coverage 771 . D increases monotonically with 771 . 

FIG. 3. BC functions for zero-width rods, (a) n{ = 0.1, di = 0.01; Upper inset: typical configuration of rods; Lower inset: 
zoo m be tween cutoffs log(ro) = —2.2, log(ri) = —1.2. (b) -qi = 0.01, cZi = 0.001; Inset: zoom with agreement on slope between 
Eq. (J3 . 2|) and the linear regression to 2 significant digits. 



FIG. 4. Apparent FD as a function of anisotropy parameter n [Eqs. (3.2), (3.9)], at two different ID coverages 



The FD decreases as the rods tend to be more parallel, and as the coverage becomes smaller. 

FIG. 5. BC functions for isotropically oriented rods of zero width and a power law distribution of lengths. The coverage here 
is r)\ — 0.1, the average rod length is (I) = 0.02, the number of rods is M — 62, and the number of possible sizes generated by 
Cantor set iterations is k m = 6. The inset shows the zoom between cutoffs log 10 (ro) = —1.89 and log 10 (n) = —0.94. 

FIG. 6. BC function for a random distribution of Cantor-rods on the interval (thick dashed-dotted line). Comparisons to 
the result for full rods of the same total length (solid line) as well as for a single perfect Cantor set with FD=D C (dashed 
line) are shown. The four regions observed for the Cantor-rods are (from left to right) of slope —1, slope — D c (which is the 
Cantor-set dimension determined by the internal structure of the rods), apparent FD (which is determined by the coverage and 
the distribution of rods), and slope —1 again. In this example there are M — 1000 Cantor rods of size d = 10 -5 • L each and 
the number of iterations in the construction of each rod is n = 5. 
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